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We develop a systematic self-consistent perturbative expansion for the self energy of Hubbard-like 
models. The interaction lines in the Feynman diagrams are dynamically screened by the charge 
fluctuations in the system. Although the formal expansion is exact — assuming that the model 
under the study is perturbative — only if diagrams to all orders are included, it is shown that for 
large-on-site-Coulomb-repulsion-C/ systems weak-coupling expansions to a few orders may already 
converge. In order to test the approximation at intermediate-to-high temperatures, we use the 
exact charge-fluctuation susceptibility from quantum Monte Carlo (QMC) simulation studies as 
input, which determines the exact screened interaction, and compare our results for the self energy 
to the QMC results. We also make comparisons with fluctuation-exchange (FLEX) approximation. 
We show that the screened interaction for the large-f/ system can be vanishingly small at a certain 
intermediate electron filling; and it is found that our approximation for the imaginary part of the 
one-particle self energy agrees well with the QMC results in the low energy scales at this particular 
filling. But, the usefulness of the approximation is hindered by the fact that it has the incorrect 
filling dependence when the filling deviates from this value. We also calculate the exact QMC 
Fermi surfaces for the two-dimensional (2-D) Hubbard model for several fillings. Our results near 
half filling show extreme violation of the concepts of the band theory; in fact, instead of growing, 
Fermi surface vanishes when doped toward the half-filled Mott-Hubbard insulator. Sufficiently away 
from half filling, noninteracting-like Fermi surfaces are recovered. These results combined with 
the Luttinger theorem might show that diagrammatic expansions for the nearly-half-filled Hubbard 
model are unlikely to be possible; however, the nonperturbative part of the solution seems to be less 
important as the filling gradually moves away from one half. Results for the 2-D one-band Hubbard 
model for several hole dopings are presented. Implications of this study for the high-temperature 
superconductors are also discussed. 

PACS numbers: 71.10.Fd, 71.18.-hy, 71.10.Hf, 74.72.-h 



I. INTRODUCTION 

The basic Hamiltonian for the simplest description of 
interacting electrons in a periodic potential of fixed lat- 
tice ions, the one-band Hubbard model, was widely in- 
vestigated after the discovery of high-temperature super- 
conductors.tJ In two dimensions, this model is widely ac- 
cepted to have an antifcrromagnetic ground state at half 
filling of its tight-binding band, although the tempera- 
ture evolution and the relation with the antiferromag- 
netism of_its insulating electron spectrum is poorly un- 
derstood. □ But, there is still ongoing debate about the 
nature of this model, in fact, there are various myster- 
ies, at band fillings close but not equal to one half. The 
model may have a superconducting ground state or be 
close to such an instability at these fillings; therefore, 
understanding the Hubbard model seems to be crucial 
for an understanding of the high-temperature supercon- 
ductors. 

There have been numerous approaches to the solution 
of the two-dimensional (2-D) Hubbard model. Exact di- 
agonalization studies are limited to very small lattices 



and mainly for this reason they are inconclusive .i Be- 
cause of statistical errors, exact quantum Monte Carlo 
(QMC) studies at band filhngs corresponding to that of 
the high-temperature superconductors (about 15% doped 
away from one half) are limited to temperatures no less 
than about a^ quarter of the electron- hopping energy 
of the model.u At such moderately high temperatures 
little information can be obtained about possible low- 
temperature instabilities such as superconductivity. But 
it should be reminded that QMC studies have led to 
a fairly good understanding of the insulating antifcrro- 
magnetic behavior of the half-filled model which corre- 
sponds to the undoped parent compounds of the high- 
temperature superconductors. 

To overcome the lattice size or temperature limitations 
of the available exact methods, approximate methods 
are needed. An important class of such approximations 
use diagrammatic formalisms. Since these approxima- 
tions are weak coupling in nature, it is important that 
the approximation sums up the physically important di- 
agrams that constitute the exact infinite perturbative ex- 
pansion. For example, self-consistent approximations like 
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fluctuation-exchange (FLEX) approximation, generalizes 
the physically important Hartree-Fock approximation by 
adding electron-hole and electron-electrqa pair scattering 
events to the bare Coulomb interaction^ 

In order to go beyond Hartrce-Fock- or FLEX-type ap- 
proximations, one has to employ additional diagrams. A 
nice formal approach to this is to renormalize, i.e., re- 
place with a perturbative expansion, the bare-interaction 
lines in the diagrams. But, this can cause an enormous 
numerical-calculational challenge, because now, the in- 
teraction lines, in principle, can depend on the momen- 
tum and frequency of the incoming and outgoing elec- 
trons, whereas for the bare Coulomb interaction, they 
don't. 

In this article, we present a new self-consistent di- 
agrammatic expansion for the electron self energy, in 
which the interaction lines are renormalized. The sys- 
tematic perturbative expansion sums up the skeleton di- 
agrams which exclude electron loops that represent a con- 
tribution to the charge fluctuations. Renormalization of 
the interaction lines is therefore achieved by the screening 
of the bare Coulomb interaction by the charge fluctua- 
tions. Although the approximation is exact, assuming 
that the model under the study is perturbative, only if 
the diagrams to all orders are included, it may converge 
rapidly if the screened interaction is considerably weaker 
than the bare Coulomb interaction. For the large-on-site- 
Coulomb-repulsion-C/ 2-D Hubbard model, we evaluate 
the exact screened interaction from the QMC data, and 
show that this is indeed the case near, but sufficiently 
doped away from, half filling. Then, by using this exact 
screened interaction, we evaluate the electron self energy 
up to the third order, and compare to the QMC results. 
We find that the imaginary part of the self energy at low 
energies almost converges to the exact QMC results at 
the third order. 

This article is organized as follows: In Section fl| we 



formulate our approximation. In Section III we present 
and discuss our results for the one-particle self energy 
along with the exact QMC and approximate FLEX cal- 
culations. We also calculate the exact QMC Fermi sur- 
faces of the 2-D Hubbard model and discuss them in Sec- 
tion |l^. And finally, we summarize our main findings in 
Section 0. 



II. METHOD 



where creates an electron with spin a at site on 
an N = Lx L square lattice with periodic boundary con- 



ditions, and Ui 



The lattice constant, a, is set 



to be 1. The Coulomb matrix element Vij =y(Ri — Rj) 
comprises the Hubbard on-site repulsion U = V{0) and 
an extended part, y(AR), for non-zero AR. For the 
calculations done in this article, we include only a first- 
neighbor hybridization t and the chemical potential fi 
in tij, and the on-site U in V^(AR). For the high- 
temperature-superconducting cuprates, jt. is estimated 
that t - 0.3-0.5 eV and U/t - 8-12MeI Wc wiU use 
U/t = 8 m all our calculations. 

In Fig. 0(a), we show the generic form of the exact 
screened interaction, Vg, expressed in terms of the bare 
Coulomb interaction, V, and the exact polarization in- 
sertion, P. In momentum-frequency space, 
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1-P{Q)V{Q)' 
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where Q — (Q, iQ). Fig. |^(b) shows the zeroth- and first- 
order diagram contributions to P, expressed in terms of 
Vg. This form of expression is useful because it allows 
one to write down a self-consistent approximation for Vg 
and P, given a one-particle Green's function. 
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FIG. 1. (a) Exact screened interaction. Vs. V is the bare 
Coulomb interaction, (b) Exact polarization insertion, P. 



We will develop our diagrammatic expansion for the 
electron self energy for the case of the 2-D one-band ex- 
tended Hubbard model. The expansion is valid in other 
dimensions as well and is easily generalizable to multi- 
band models. The Hamiltonian is 
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These two simply-related quantities are also similarly 
related to a physical quantity easily measurable in QMC 
studies, the charge susceptibility, Xc- The expression 
for this two-particle correlation function as a thermody- 
namic expectation value is 

x(An,(T)An,(0)), (3) 
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where ni(r) = J2cr 4<ji'^)^i'^('^)^ ^"^'^ ^ni{T) = mir) - 
K(r)). 

In principle, the definition of a polarization insertion, 
P, is arbitrary, and so is of a screened interaction, V^, 
through Eq. (0) . P defined in this article through the di- 
agrammatic expansion in Fig. |l|(b) actually corresponds 
to the polarization insertion of the charge susceptibility, 
Xc- The exact Xc is expressible as a geometric series in 
terms of the exact P, such that 

" ^'^^^ " 1 - PiQMQ) 

(see also Fig. ^). Combining Eqs. and (^), we obtain 
an important expression for the screened interaction, Vg, 
in terms of the physical charge susceptibility of the sys- 
tem, Xc; 

V,{Q) = V{Q) - V{Q)xciQ)ViQ). (5) 



FIG. 2. The charge susceptibility, Xc, expanded in terms of 
the polarization insertion, P. V is the bare Coulomb interac- 
tion. 

At this point, we are ready to calculate the exact Vg by 
using the exact Xc obtained by QMC. Using V{Q) = U 
for the Hubbard model, we have, 

Vg[Q) = U -U\,{Q). (6) 

One can also Fourier transform Vs{Q) to the real space, 
by using, 

V^s(AR, ^^) = ]^Y. ^''^•^''^(Q, i^). (7) 
Q 

Again, for the Hubbard model, 

V;(AR, m) = C/(5(AR) - [/2xc(AR, iil), (8) 

where S is the Kronecker delta. In Fig. ||, we plot 
the on-site and time-independent component of the ex- 
act screened interaction obtained from QMC, T4(AR = 
O^ifl = 0), for different hole dopings, a; = 1 — (n), mea- 
sured with respect to half filling. Although the bare (un- 
screened) Coulomb interaction, U, is equal to 8t, strik- 
ingly, the corresponding component of the exact screened 
interaction is substantially weaker for intermediate hole 
dopings, X, actually vanishing and changing sign, i.e., be- 
coming attractive, ft a; ~ 15%. We will refer to this effect 
as "overscreening,"l3 and explain below why it happens. 

From Eq. (||), because of the term, one can sus- 
pect that for sufficiently large U, this term will win and 
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FIG. 3. V;(AR = 0,in = 0) as a function of the hole 
doping, X = 1 — (n), obtained from QMC for the 2-D Hubbard 
model. Parameters are U/t = S\ T/t = 1/3, where T is the 
temperature; and the lattice size is 8 x 8. 
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V^(AR = 0, if] = 0) will become negative (attractive), 
since Xc(AR = 0,ifl = 0) is always positive and ap- 
proaches a nonzero JJ-independent limit [of 0{x/t) by 
a simple U = oo scaling argument] for U ^ t. On 
the other hand, at half filling, charge fluctuations are 
suppressed; so that Xc(AR, zfi = 0) ~ 0{t'^/U^) and 
14(AR = 0,ir2 = 0) ~ U, for U :^ t. Hence, presence 
of a sufficiently large U and finite hole doping off half 
filling, X, is necessary for overscreening. Also, note that 
overscreening is not restricted to the on-site component 
of the Coulomb interaction; the nearest-neighbor compo- 
nents, and so forth, can also become attractive because 
of the overscreening effect HI 

At this point, we are ready to write down a controlled 
perturbative expansion for the one-particle self energy. 
Since the exact screened interaction, Vs, was shown to be 
substantially weaker than the bare interaction, the Hub- 
bard U, one can expect that an expansion in terms of Vs 
should converge much more rapidly than a brute-force 
expansion in terms of U. As a matter of fact, it is not 
difficult to pick up the diagrams which contribute to this 
expansion. Keeping in mind that our approximation for 
the self energy will be self consistent in the sense that all 
one-particle-Green's-function lines will correspond to the 
full Green's function, we should allow only the skeleton 
diagrams in which none of these lines have any self-energy 
diagrams explicitly attached to them. As a second step, 
bearing in mind the fact that the interaction lines in our 
diagrams, which correspond to Vg, already include the 
polarization diagrams in them [Fig. 0(a)],' we omit any 
self-energy diagram in which an interaction line has a 
polarization diagram [see Fig. 0(b)] inserted. These are 
basically the only rules needed in the diagram selection. 
If all the diagrams to infinite order are included, an exact 
perturbation series will result, like in any other rigorous 
perturbative expansion. Note that the interaction lines in 
our diagrams, or V^, represent the bare Coulomb interac- 
tion, U, screened by the charge fluctuations, Xc [Eq- @]- 
Therefore, we have derived a diagrammatic expansion in 
terms of the charge fluctuations, which is exact if carried 
out to all orders. Moreover, because of the overscreening 
effect discovered by the analysis of the interaction lines 
(Vs) obtained by QMC exactly, a weak-coupling expan- 
sion to a first few orders is expected to converge in a 
rapid, controlled fashion. 

In Fig. 0, we show all the self-energy diagrams up to 
the third order. The Hartree diagram in Fig. 0(a) is 
written separately in terms of the bare interaction (U in 
this case) to prevent double counting; for the one-band 
Hubbard model, it is a trivial constant which is equal to 
U{n). The other first order diagram in Fig. 0(a) is like 
the Fock diagram with the bare interaction line replaced 
by the screened interaction. This diagram combines the 
bare-interaction Fock diagram (which is again equal to a 
trivial constant, —U{n)/2, in this case) with the charge- 
fluctuation self energy diagram found in FLEX-like ap- 
proximations; because, the screened interaction is the 
sum of the bare interaction and the charge-fluctuation 



propagator [see, e.g., Eq. (0) or (0)]. In the second order, 
there is only one diagram contributing to the self energy 
[Fig. 0(b)]; and in the third order, the number of all the 
contributing diagrams is six [Fig. 0(c)]. 
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FIG. 4. All the contributing self energy diagrams in 
the first (a), second (b), and third order (c) of the 
screened-interaction (Vs) expansion. The weak-coupling ex- 
pansion would become exact (assuming the model under 
study is perturbative) when carried out to all orders. 

At this point, we would like to discuss the fact that, 
the screened-interaction expansion developed here is a 
low-energy approximation. In Fig. 0, we plot the exact 
screened interaction, V^(Q,i51), which was obtained by 
QMC, as a function of the frequency, il. For the chosen 
Q point, Vs is actually attractive and small at zero fre- 
quency, and vanishes at a slightly higher frequency. But 
as we approach the first Matsubara frequency, Vs rapidly 
grows, and it is already more than half of the bare, un- 
screened Coulomb repulsion. As the frequency increases 
further, Vs finally approaches to the bare repulsive inter- 
action, U = 8t. It is now clear that the high-energy part 
of the Vs is not weak, and our approximation is limited 
to the low-energy scales. Therefore, for the tempera- 
ture studied here (T — t/3), the calculated self energy 
is expected to be accurate only for the first Matsubara 
frequency, ttT k t, a small-enough fraction of the width 
of the tight-binding band, which governs the energy scale 
for Vs. 
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FIG. 5. The exact screened interaction, Vs(Q,if2), calcu- 
lated by QMC, shown as a function of the frequency, Q.. The 
momentum point is Q = (7r/2, 0) on an 8 x 8 lattice. Param- 
eters are U/t = 8; T/t — 1/3, where T is the temperature; 
and X = 14.0%, where x = 1 — (n). 



= ViQ = 0){n)--J2 ^(Q' k)V,{Q)G{k - Q), 



(10) 



where G{k) is the one-particle Green's function and k — 
(k, zcj). Eq. ^ immediately implies that for a given Q, 
if P{Q) = oo, then A(Q, k) — oo, because G(/c) is always 
finite. Even though A((3, k) in Eq. ( [lO| ) may diverge, the 
one-particle self energy, which is the actual physical 

quantity, remains finite; because, if A{Q, k) diverges at 
a given Q, then Vs{Q) will vanish, therefore canceling 
the divergence. But, there is still a concern that a weak- 
coupling expansion for A in Fig. |^(b), which generates the 
diagrams in Fig. ^, may not be sufficient because of this 
diverging behavior. We will put our theory to numerical 
testing in the next section. 
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Before finishing this section, we would like to discuss 
a peculiar property of the irreducible polarization inser- 
tion and the associated three-point vertex function. As 
we discussed earlier, the exact screened interaction, T4, 
as obtained from QMC, vanishes at particular (AR, iVl) 
points, and it will vanish at particular (Q,ir2) points 
likewise; this was the motivation behind developing a 
weak-coupling expansion in terms of Vs{Q). But from 
Eq. (^), Vs{Q) = immediately imphes that P{Q) = oo 
at the same Q — {Cl,iQ) point. This divergence in the 
polarization insertion of the charge susceptibility does 
not imply a divergence in the charge susceptibility it- 
self, which is the actual physical quantity. In fact, from 
Eq. (|), P{Q) = oo simply implies that Xc{Q) = l/V{Q). 
Although this divergence is not a worry from a physics 
point of view, it imposes technical difficulties on the ex- 
pansion of P{Q) [see Fig. |](b)] along these singularities; 
to overcome this difficulty in an actual calculation of 
P{Q), it would be more appropriate to expand its inverse, 
1/P{Q), rather than to expand P{Q) itself directly. 

Both the polarization insertion and the self energy can 
be expressed in terms of a three-point irreducible vertex 
function, A{Q, fc), by the exact relations (see Fig. ||), 
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FIG. 6. The polarization insertion, P(Q) (a), and the 
one-particle self energy, E(A;) (b), in terms of the three-point 
irreducible vertex, K(Q,k). G{k) is the one-particle Green's 
function, Q — (Q,if2), and k — (k,itj). V and 14 are the bare 
and the screened interaction respectively. 



III. ONE-PARTICLE SELF ENERGY 

We perform our numerical calculations for the self en- 
ergy on periodic discreet lattices of sizes 4x4 and 8x8, 
so that they can be compared with the QMC results. We 
use the Matsubara-frequency representation, and employ 
a fermion frequency cutoff of ^ 35t. All the QMC calcu- 
lations were done at a temperature, T = i/3, and with 
the cutoff employed here, this corresponds to 16 positive 
Matsubara frequencies. Note that QMC calculations are 
performed in the imaginary time, and a time discretiza- 
tion is made for numerical purposes. The number of time 
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slices in [0,/3 = l/T] used in the QMC calculations was 
80, which is equivalent to using 40 positive Matsubara 
frequencies — even larger a cutoff than our diagrammatic 
calculations. We chose a smaller cutoff for our diagram- 
matic calculations to avoid the high-frequency errors in 
the Vs obtained from QMC, which goes directly into the 
calculations. We vary the chemical potential, /i, in or- 
der to match the filling, (n), with that of the QMC. In 
order to calculate (n), the one-particle Green's function, 
G'(k, iuj), has to be summed over the Matsubara frequen- 
cies, and to prevent cutoff effects in this relatively-slow- 
converging sum, we effectively extend the sum to infinity 
by using the standard — replace the high-frequency part 
of the sum with the result for the non-interacting case — 
trick. 

In Fig. 1^, we plot the real and imaginary parts of the 
self energy as a function of Matsubara frequencies along 
with the results from QMC calculations. We chose a 
particular momentum point close to the Fermi surface, 
k = (tt, 0), but the k dependence of both the QMC and 
the screened-interaction-expansion calculations near the 
Fermi surface is small. QMC error bars are omitted in all 
but one plot in this article for clarity, but the estimated 
QMC error bar for the real and imaginary parts of the self 
energy at the first frequency is about zLO.lt. The corre- 
sponding error bar for the QMC Green's function is about 
±0.02t^^. Also, the error bars for the self energy grow 
with the frequency, whereas, they are roughly constant 
for the Green's function. Note that, the temperature in 
these calculations is T — t/3, and the second Matsubara 
Frequency, SttT, is already about 5.2t, which is a substan- 
tial fraction of the (non-interacting) bandwidth of our 
model, 8t. For the reasons we explained in the previous 
section, our approximation is a low-energy approxima- 
tion, which is unable to capture this regime. But, the first 
Matsubara frequency, ttT « l.Ot, is low enough compared 
to the bandwidth, and remarkably, at the 14% hole dop- 
ing, where the screened interaction almost vanishes (see 
Fig. 1^), our approximation for the imaginary part of the 
self energy almost converges to the QMC result within 
the error bar [Fig. 0(b)]. This shows that our expan- 
sion may indeed be converging vary rapidly at the hole 
doping at which the on-site component of the screened 
interaction at zero frequency, T4(AR = 0,ifl = 0), van- 
ishes. But, when the hole doping is varied away from 
this point, Im S varies in opposite direction with respect 
to the exact QMC results. Naturally, having the correct 
doping dependence is one of the most desired aspect of 
any approximation, which, therefore, substantially limits 
the usefulness of this expansion. It is quite possible that 
the reason for this failure may be intimately related to 
the vertex divergence problem discussed in the previous 
section. In Fig. |[ we show the evolution of the self energy 
with the order of the diagrams included. The imaginary 
part, in fact, seems to be converging to the QMC result 
at the lowest frequency with the increasing order, and it 
is about the same as the QMC result within the QMC 
error bar (see Fig. |l^) at the third order. But the real 



1.5 



zi. 

I 



-0.5 



-0.5 



-1 



-1.5 



-2 



-2.5 



(a) 
















▲ 


,0, 




— ■ * 


— • 


• 







" 


■ 






A' / 
0' /' 


n' ' 



4 6 
CO/ t 



10 



(b) 

























10 



CO/ t 



FIG. 7. The one-particle self energy as a function of the 
Matsubara frequency, shown for the screened-interaction ex- 
pansion (solid lines and filled symbols) along with exact (ex- 
cept for statistical errors) QMC results (dashed lines with 
empty symbols). The k point is (7r,0). Hole dopings, 
X = \ — (n), are 10.0 (squares), 14.2 (circles), and 17.5% 
(triangles). Other parameters are U jt = % and T jt — 1/3. 
Lattice sizes are 4x4 for the screened-interaction expansion 
and 8x8 for QMC. Hole dopings for both the expansion and 
the QMC correspond to averages on 4 x 4 k-space meshes. 
Actual QMC hole dopings averaged on 8 x 8 k-space meshes 
are 9.7, 14.0, and 17.1%. The real (a) and the imaginary (b) 
part. [We use 8x8 QMC results in these comparisons because 
QMC calculations were performed only on 8 x 8 lattices for 
most hole dopings; but the difference between 4x4 and 8x8 
QMC results are well within the error bars (see Fig. ^ and 
the associated discussion at the end of this section).] 
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part seems to be going in the wrong direction with the 
increasing order. 

The real part of the self energy which was calculated 
by the screened-interaction expansion actually does not 
do very well against the exact QMC results [Fig. 0(a)] 
at all, but, as we will analyze below, this is not a failure 
of the expansion itself. In fact, this is merely the quan- 
tity which determines the location of the Fermi surface. 
Luttinger theoremEi states that, for any diagrammatic ex- 
pansion, the total volume (area in 2-D) enclosed by the 
Fermi surface is equal to the total filling, (n); therefore 
our approximation will obey this statement. But then, 
this simply implies that the exact QMC results may not 
obey the Luttinger theorem for the volume of the Fermi 
surface. To investigate this further, we calculated the ex- 
act QMC Fermi surfaces for various dopings, which will 
be discussed later. 

Fig. ^ shows results for the FLEX calculations in com- 
parison with the QMC. FLEX approximation employed 
here uses exchange of charge and spin fluctuations, i.e., 
the particle-hole-channel, but omits the particle-particle 
fluctuations (for the Hubbard model there would only 
be the particle-particle fluctuations of the singlet type). 
The results for the real part are very similar to the 
screened-interaction expansion. FLEX has a somewhat 
better qualitative behavior at high frequencies. Low fre- 
quency behavior of the real part for FLEX and screened- 
interaction expansions are expected to be very similar be- 
cause they are both expected to obey the Luttinger the- 
orem. The real part for FLEX at the lowest frequency is 
doing somewhat worse than the screened-interaction ex- 
pansion, probably because Im S for FLEX is somewhat 
bigger, indicating that FLEX results are more distant to 
their zero- or low-temperature values [compare Figs. 0(b) 
and ^(b)]. Im S is expected to vanish at zero tempera- 
ture and zero frequency also according to the Luttinger 
theorem. FLEX results for the Im E show a correct qual- 
itative behavior at high frequencies. Although they also 
capture the doping dependence correctly, the magnitude 
of the doping dependence is smaller by an order. In addi- 
tion, they don't do well quantitatively at any particular 
doping. 

The screened-interaction-expansion calculations per- 
formed in this article are computationally very time con- 
suming. The reason for this is the requirement for the 
summation over many frequency and momentum vari- 
ables for the third and higher order diagrams. In order to 
go to the fourth order, one would probably need Monte 
Carlo summation techniques. Even at the third order, 
we were able to perform our calculations on 4 x 4 lat- 
tices only. In Fig. 10 we make some comparisons of the 
calculations on different lattice sizes. 4x4 and 8x8 
QMC results are compared to the screened-interaction- 
expansion results. FLEX results are shown as well. The 
difference between the results for the two lattice sizes 
seems to be very small. Estimated QMC error bars for 
the Im E on the 8x8 lattice are also shown. These are 
the approximate upper bounds to the statistical errors 
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FIG. 8. The one-particle self energy as a function of the 
Matsubara frequency, shown for the 1st- (triangles), 2nd- 
(squares), and 3rd-order (circles) screened-interaction expan- 
sion (solid lines and filled symbols) along with exact (except 
for statistical errors) QMC results (dashed lines and empty 
circles). The k point is (7r,0). Hole doping, x = 1 ~ (n), is 
14.2%. All parameters are the same as in Fig. ^. The real (a) 
and the imaginary (b) part. 
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associated with the QMC calculations. 
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FIG. 9. The one-particle self energy as a function of the 
Matsubara frequency, shown for FLEX approximation (solid 
lines and filled symbols) along with exact (except for statisti- 
cal errors) QMC results (dashed lines with empty symbols). 
The k point is (7r,0). Hole dopings, i = 1 — (n), are 10.0 
(squares), 14.2 (circles), and 17.5% (triangles). Other param- 
eters are U/t = 8 and T/t = 1/3. Lattice sizes are 4x4 for 
FLEX and 8 X 8 for QMC. Hole dopings for both the FLEX 
and the QMC correspond to averages on 4 x 4 k-space meshes. 
Actual QMC hole dopings averaged on 8 x 8 k-space meshes 
are 9.7, 14.0, and 17.1%. The real (a) and the imaginary (b) 
part. 



IV. EXACT QMC FERMI SURFACE 

In this section, we will discuss the exact QMC Fermi 
surfaces calculated. Although there are several equiva- 
lent definitions for the Fermi surface, the definition which 
locates the Fermi surface as the zero frequency poles of 
the Green's function is most suitable for finite tempera- 
tures. One finds the solution of e(k)+Re S(k, uj — Q) — jj, 
in order to locate such poles.il For a Fermi liquid, these 
poles correspond to the quasiparticles at the Fermi sur- 
face. But, the definition is general, and can also be used 
at finite temperatures. Note that, the Luttinger theorem, 
which relates the volume enclosed by the Fermi surface 
to the total electron filling, is strictly valid only at zero 
temperature; but, as long as the temperature is not un- 
reasonably high, one finds that the Luttinger theorem is 
satisfied in an approximate manner, and the only change 
to the quasiparticlc picture is some temperature broaden- 
ing due to the finite Im I](k, w = 0). As mentioned in the 
previous section, the quantity, ImE(k, w = 0), is indeed 
a good indication of how close one is to T = 0, which 
should actually vanish at T = for a system obeying the 
Luttinger theorem. In order to numerically calculate the 
Fermi surface on a finite 8x8 QMC lattice, we first find 
the solution of £(k)-|-Re S(k, = 0) = /i by interpolating 
between the k points. Moreover, S(k, a; — 0), which is 
very close in value to the one at the lowest available fre- 
quency, I](k, w = iT^T), is found by linear extrapolation 
using the lowest-two frequencies. Due to this numerical 
approximation procedures, the systematic interpolation 
and extrapolation errors up to a few percent are possible. 

Fig. |ri](a) shows the exact QMC Fermi surfaces for 
actual hole dopings, x = 1 — (n), ranging from 5.0 to 
33.5%. These actual hole dopings in decreasing order are 
X = 33.5 [smallest surface centered around k = (0,0)], 
23.5, 20.3, 17.1, 14.0, 9.7, and 5.0% (the arcs centered 
around the Brillouin-zone corners). We found the ar- 
eas enclosed by the Fermi surfaces by numerical integra- 
tion, and calculated the hole dopings implied by the Lut- 
tinger theorem corresponding to these areas. Note that, 
due to the curve-fitting procedures employed for the dis- 
crete (k, Lo) space, these areas are uncertain within a few 
percent. The results in the same order as the actual 
hole dopings are 33.5, 20.5, 14.1, 7.1, -4.3, -19.8, and 
—32.8%. We plot these values against the actual dop- 
ing in Fig. 0(b). Considering the fact that we are at 
a moderately high temperature, T — t/3, and account- 
ing for the curve-fitting errors, the Luttinger theorem 
seems to be satisfied at least qualitatively for hole dop- 
ings, X > 17%, i.e., the so-called overdoped region in 
the context of high-temperature superconductivity. On 
the other hand, in the underdoped region, QMC Fermi 
surfaces start deviating from the Luttinger theorem in a 
systematic qualitative matter. As one proceeds towards 
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FIG. 10. The one-particle self energy plotted as a func- 
tion of the Matsubara frequency for the purpose of eval- 
uating lattice-size effects. Results for QMC (squares), 
screened- interaction expansion (circles), and FLEX approx- 
imation (triangles) are shown. The k point is (7r,0). Sohd 
lines and filled symbols represent 8x8 lattices for QMC and 
4x4 otherwise. Dashed lines and empty symbols represent 
4x4 lattices for all. Parameters are C//t = 8 and T jt — 1/3. 
Hole dopings, x = 1 — (n), averaged on 4 x 4 k-space meshes 
are 14.2% (solid lines and filled symbols) and 13.8% (dashed 
lines and empty symbols). Actual, 8 x 8- mesh averaged hole 
doping for the 8 x 8 QMC lattice is 14.0% (sohd line with filled 
squares). Note that the 4 x 4-mesh averaged hole-doping val- 
ues for all data sets match each other among the solid and 
the dashed lines. The real (a) and the imaginary (b) part. 
Also shown are the QMC error bars for the Im E on the 8x8 
lattice. 



half filling, Fermi surface starts shrinking towards the 
k = (7r,7r) point, looking like an electron-doped system. 
This is completely inconsistent with the concepts of the 
band theory and Fermi liquid. Fermi surface of the Hub- 
bard model gradually gets smaller and disappears when 
the hole doping gets smaller in the underdoped region. 
Therefore, this seems to be a doping induced transition 
between a Fermi-liquid metal and a strongly-correlated 
insulator. This extreme violation of the Luttinger theo- 
rem may imply that there is no asymptotically conver- 
gent diagrammatic expansion for the Hubbard model in 
the underdoped region. But, it looks like, as the hole 
doping increases past the optimal doping 15%), the 
system looks more like a conventional Fermi liquid, and 
diagrammatic expansions seem more likely to be possible. 
In Fig. |ll](b), we plot the hole dopings, x, deduced from 
the area enclosed by the Fermi surface according to the 
Luttinger theorem against the actual hole dopings. Sys- 
tematic deviation towards half filling is clearly seen. Note 
that the line is only a guide to the eye; it is quite possible 
that the area of the QMC Fermi surface at dopings very 
close to half filling gets vanishingly small. We don't have 
data for such dopings at the point, so, further investiga- 
tion is neeassary to determine the existence of such "hole 
pockets. "Ily In order to rule out temperature effects one 
might also study lower temperatures, which could be ac- 
cessible for such dopings very close to half filling. This 
figure also shows the results for FLEX calculations at the 
same temperature. Since FLEX is expected to satisfy the 
Luttinger theorem at low temperatures, the deviations in 
this case are known to be mainly because of temperature 
effects. But, these deviations are much more reasonable, 
and the FLEX results show approximately correct slope 
for all dopings. Doping deduced from the Luttinger the- 
orem for FLEX calculations also have the correct sign 
for all dopings. We show the FLEX Fermi surfaces at 
5.0% doping for the temperature studied above (T = t/3) 
along with a low-temperature Fermi surface in Fig. n2. 
The low- and intermediate-temperature Fermi surfaces 
for FLEX look very similar, and the low-temperature 
case is almost identical to the noninteracting case (for the 
Hubbard model, which is also the same as the Hartree- 
Fock case). The comprehensive results we obtained for 
the Fermi surface are also qualitatively consistent with 
the other recent QMC study of the JIubbardD and high- 
temperature expansion study of t-Jt3 models. 



V. SUMMARY AND CONCLUSIONS 

In this article, we studied a new approximation for the 
one-particle self energy of the large-?/ Hubbard-like mod- 
els. An expansion was written out in terms of a screened 
interaction, which is nonlocal and retarded, although it 
acts between local and instantaneous charge densities. 
The screening is achieved by the charge fluctuations. The 
expansion is an exact perturbative expansion, although 
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FIG. 11. (a) The exact QMC Fermi surfaces for various 
hole dopings. Parameters are U/t = 8 and T/t = 1/3. 
Actual hole dopings, x = 1 — (n), in decreasing order are 
X = 33.5 [smallest surface centered around k = (0,0)], 23.5, 
20.3, 17.1, 14.0, 9.7, and 5.0% (the arcs centered around the 
Brillouin-zonc corners). The lattice size is 8 x 8. (b) The hole 
dopings, X, deduced from the areas enclosed by the Fermi sur- 
faces according to the Luttinger theorem, plotted against the 
actual hole dopings, x = 1 — {n) (circles and the solid line). 
Also shown are the FLEX results at the same temperature 
(squares and the dashed line). The long-dashed line without 
symbols corresponds to the Luttinger theorem. The lattice 
size for the FLEX calculations is 16 x 16. 




-0.5 



FIG. 12. The FLEX Fermi surface at the temperature 
T/t =1/3 (long-dashed line) and at a lower temperature, 
T/t = 1/16 (solid line). This is the lowest available tem- 
perature for the lattice size studied (16 x 16) at this doping. 
Parameters arc U/t — 8 and x = 1 — (n) = 5.0%. Hole dop- 
ings, X, deduced from the areas according to the Luttinger 
theorem are 7.6% {T/t = 1/16) and 9.5% {T/t = 1/3). Note 
that the low-temperature FLEX Fermi surface is almost the 
same as the noninteracting (or Hartree-Fock) Fermi surface 
(not shown). 
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this doesn't necessarily mean that it is an exact solu- 
tion because of the possibihty of nonperturbative nature 
of the large- 1/ systems. We carried out the expansion 
up to the third order. The screened interaction was de- 
termined by QMC calculations exactly, and then directly 
put into the expansion calculations. The local component 
of the exact screened interaction, which screens the bare 
Hubbard C/, was found to be vanishing at certain small 
frequencies at hole dopings near 15%. This is a very en- 
couraging behavior for a weak-coupling expansion, which 
indicates that such expansions might be possible at low- 
enough frequencies where the screened-interaction is very 
small. However, the analysis of the three-point-vertex 
function, which, together with the screened interaction, 
determines the one-particle self energy, showed that it 
diverges at the same frequency-momentum points where 
the screened interaction vanishes. The expansion for the 
self energy had only limited success against the exact 
QMC results. Although the expansion gave excellent re- 
sults for the low-energy part of Im E near 15% hole dop- 
ing, where the screened interaction vanishes, it had the 
wrong doping dependence around this same value. More- 
over, the real part of the QMC results were not consis- 
tent with a diagrammatic expansion which should obey 
the Luttinger theorem, although, we can't rule out the 
possibility of temperature effects. If it were because of 
temperature effects, this would also mean that our ap- 
proximation doesn't work for Re E at all, since with the 
increasing order, results for the expansion moved in the 
wrong direction. In contrast, Im E showed a converg- 
ing behavior near 15% doping, as mentioned above. We 
also calculated FLEX self energies, which had the cor- 
rect qualitative doping dependence and high-frequency 
behavior, but did actually worse quantitatively near 15% 
doping. 

In order to investigate the QMC Re E results, we nu- 
merically calculated the exact QMC Fermi surfaces. The 
results showed extreme violation of the Luttinger theo- 
rem near half filling, which might be because of a doping- 
induced metal-insulator transition arising from strong- 
coupling effects. The Fermi surface of the Hubbard model 
seemed to be reducing to hole pockets centered around 
k = (tt, tt) at small hole dopings near half filling. Fur- 
ther QMC investigation at smaller hole dopings and lower 
temperatures might bring more evidence into this. 

As far as the high-temperature superconductors are 
concerned, exact QMC results might be indicating that 
the underdoped region of these materials is a transition 
between a strong-coupling Mott-Hubbard insulator and 
a Fermi metal. The disappearance of the Fermi surface 
with the decreasing hole doping can also be viewed as 
an opening of a pseudogap.EJ For the full investigation 
of this region at low temperatures, other nonperturba- 
tive treatments might be necessary. However, the partial 
success of the screened-interaction expansion developed 
here shows that it might be useful near optimal doping, as 
well as other diagrammatic expansions in the optimally- 
doped and overdoped regions. 



ACKNOWLEDGMENTS 

G. E. acknowledges useful discussions with Richard T. 
Scalettar and Andrew K. McMahan. This work was 
partially supported by National Science Foundation un- 
der the Grants No. DMR-92-15123 and No. DMR-99- 
70291. Work at University of California, Davis was sup- 
ported in part by an Accelerated Strategic Computing 
Initiative Grant and by Materials Research Institute of 
Lawrence Livermore National Laboratory. Computing 
support from University Computing and Networking Ser- 
vices at University of Georgia is gratefully acknowledged. 



* Present address: Department of Physics, University of Cal- 
ifornia, Davis, California 95616-8677; and Lawrence Liver- 
more National Laboratory, University of California, Liver- 
more, California 94550-9234 

^ E-mail address: esirgen@ucdavis.edu 

* Present address: Institut fiir Theoretische Physik, Technis- 
che Universitat Graz, 8010 Graz, Austria 

^ E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

^ M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 
^ W. Fettes, I. Morgenstern, and T. Husslein, Computer 

Physics Communications 106, 1 (1997); H. Q. Lin, J. E. 

Hirsch, and D. J. Scalapino, Phys. Rev. B 37, 7359 (1988). 

C. Grober, R. Eder, and W. Hanke, Phys. Rev. B 62, 4336 
(2000). 

^ G. Esirgen, H.-B. Schiittler, and N. E. Bickers, Phys. Rev. 

Lett. 82, 1217 (1999); and references therein. 
^ M. S. Hybertsen, E. B. Stechel, M. Schliiter, and D. R. 

Jennison, Phys. Rev. B 41, 11068 (1990); C.-X. Chen and 

H.-B. Schiittler, thid. 43, 3771 (1991); S. B. Bacci, E. R. 

Gaghano, R. M. Martin, and J. F. Annett, ihid. 44, 7504 

(1991) . 

H.-B. Schiittler and A. J. Pedro, Phys. Rev. B 45, 7588 

(1992) . 

^ H.-B. Schiittler, C . Grober, H. G. Evertz, and W. Hanke, 
cond-mat/9805133| ; 0104300 (to be published). 

^ J. M. Luttinger, Phys. Rev. 119, 1153 (1960). 
^° A. P. Kampf and J. R. Schrieffer, Phys. Rev. B 42, 7967 
(1990); S. A. Trugman, Phys. Rev. Lett. 65, 500 (1990); 

D. Duffy and A. Moreo, Phys. Rev. B 51, 11882 (1995). 
" W. O. Putikka, M. U. Luchini, R. R. P. Singh, Phys. Rev. 

Lett. 81, 2966 (1998). 

For experimental evidence of the pseudogap, see, e.g., the 
following and references therein: Z.-X. Shen and J. R. 
Schrieffer, Phys. Rev. Lett. 78, 1771 (1997) [ARPES]; R. 
Nemetschek, M. Opel, C. Hoffmann, P. F. Miiller, R. Hackl, 
H. Berger, L. Forro, A. Erb, and E. Walker, Phys. Rev. 
Lett. 78, 4837 (1997) [Raman]; G. V. M. Williams, J. L. 
Tallon, E. M. Haines, R. Michalak, and R. Dupree, Phys. 
Rev. Lett. 78, 721 (1997) [NMR]. 



11 



